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A computationally tractable version of the Bohr-Mottelson collective model is presented which 
makes it possible to diagonalize realistic collective models and obtain convergent results in relatively 
small appropriately chosen subspaces of the collective model Hilbert space. Special features of the 
proposed model is that it makes use of the beta wave functions given analytically by the softened-beta 
version of the Wilets-Jean model, proposed by Elliott et al., and a simple algorithm for computing 
SO(5) D SO(3) spherical harmonics. The latter has much in common with the methods of Chacon, 
Moshinsky, and Sharp but is conceptually and computationally simpler. Results are presented 
for collective models ranging from the sherical vibrator to the Wilets-Jean and axially symmetric 
rotor-vibrator models. 
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D , This paper presents a computationally tractable version of the standard Bohr-Mottelson (BM) collective model 
Q ' pi 0, • The need for such a version of the model arises because the expansion of rotational wave functions in a 
. , spherical vibrational basis is so slowly convergent that the diagonalization of a general coupled- rotor-vibrator collective 
■ model Hamiltonian in this basis is virtually impossible. The proposed method constructs a basis for the collective 
model in which beta wave functions, centred about a non-zero equilibrium value, are given analytically following the 
H ■ methods of Elliott et al. Q and ref. [j| , and new methods are developed for computing the complementary rotational 
' and gamma wave functions. The motivation for this development is to be seen in the context of a sequence of steps from 
ON , a phenomenological model description of nuclear collective structure to a microscopic interpretation of what it means. 
I follow the strategy of giving a phenomenological collective model a microscopic foundation by first formulating it in 
algebraic terms and subsequently constructing representations of its algebraic expression on a microscopic shell model 
Hilbert space. 

[ Early advances in the shell model theory of collectivity came with the identification of seniority and symplectic 
symmetry with pairing 0,0,0 and su(3) with rotations [jj. An early description of the use of symmetry in nuclear 
""^ ' structure was provided by Parikh [Toj . 

Algebraic methods have subsequently had an enormous influence on collective model theory (cf. [Tl| for a review). 
1 1 ' Two influential developments have been: the Frankfurt version of the collective model and the Interacting 

Boson Model (IBM) [lj|. The Frankfurt methods were developed to give solutions to the collective model in the 



intermediate region between the analytically solvable vibrational and rotational limits. They made major use of the 
\ algebraic structures associated with the five-dimensional harmonic oscillator. The IBM achieved a reduction of the 
Frankfurt program to a finite-dimensional space by compactifying the algebraic structure of the collective model to 
the U(6) symmetry group of the six-dimensional harmonic oscillator and restricting consideration to single (finite- 
dimensional) U(6) irreps. The IBM has three exactly solvable limits, corresponding to similar limits of the BM 
' model. 

Major developments towards the goal of formulating a microscopic (shell model) theory of collective states have 
also followed the algebraic approach (a historical survey was provided in the review article of ref. EH). The latter 
developments were based on the symplectic model |l6j . 

The symplectic model and its several variations have not been applied widely to fit detailed nuclear data because 
that was not their purpose (cf. refs. [H[l3 for some recent applications). Their purpose was to obtain a fundamental 
explanation of nuclear collective dynamics and the way it emerges from interacting nucleons. Thus, even if one were 
to succeed in formulating a completely satisfactory theory of nuclear collective dynamics, it would hardly be a simple 
theory. Nor would it serve the purpose of every-day analysis of nuclear data. For this purpose, one continues to need 
simple phenomenological models, albeit preferably ones with a microscopic foundation. With this concern in mind, 
we revisit the BM collective model with the benefit of insights acquired from the Frankfurt model, the IBM, and the 
symplectic model. 



II. THE STANDARD SOLVABLE LIMITS OF THE COLLECTIVE MODEL 



The classic BM collective model shares a Hilbert space with the five-dimensional harmonic oscillator. Indeed, 
in its harmonic vibrational limit, its spectrum and eigenstates are precisely those of the five-dimensional harmonic 
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oscillator. Thus, it has rich algebraic, geometrical, and analytical structures all which arc exploited in this paper. 
It has a spectrum generating algebra given by the semi-direct sum Lie algebra [HW(5)]u(5) and a corresponding 
dynamical group [HW(5)]U(5). The Lie algebra [hw(5)]u(5) is spanned by five (L = 2) pairs of d-boson (phonon) 
operators {d^ , dl; v = 0,±1,±2} and the infinitesimal generators {d^d v } of U(5), where the d-bosons satisfy the 
commutation relations [cPjdJ,] = 8^. 

The collective model has three well-known algebraically solvable limits: the harmonic vibrator model, the Wilcts- 
Jean (gamma- soft) model [3, and the axially-symmetric rigid-rotor model. These submodels are associated with 
dynamical subgroup chains corresponding to different paths through the set of groups 

pW(5)]U(5) D [R 5 ]SO(5) D [R 5 ]SO(3) 

U U U (1) 

U(5) D SO(5) D SO(3) 

starting with [HW(5)]U(5) and ending with SO(3), where R 5 is the group with Lie algebra spanned by the quadrupolc 
moments 

[Q„ = ±={dl + d v ),v = Q,±l,±2] (2) 

with d v = {-l) v d~ v . Thus, 

[HW(5)]U(5) D U(5) D SO(5) D SO(3) (3) 

is a dynamical subgroup chain for the harmonic vibrator model, 

[HW(5)]U(5) D [R 5 ]SO(5) D SO(5) D SO(3) (4) 

is a dynamical subgroup chain for the Wilets-Jean (beta-rigid, gamma-soft) model, and 

[HW(5)]U(5) D [R 5 ]SO(5) D [R 5 ]SO(3) D SO(3) (5) 

is a dynamical chain for the rigid-rotor (beta- and gamma-rigid) model. The above subgroup chains are discussed in 
more detail, for example, in ref. |TT1 |. 

The solution of more general collective model Hamiltonians has been tackled, for example, by Hess et al. 01 m 
the U(5) D SO (5) D SO (3) basis for the collective model Hubert space of Chacon, Moshinsky and Sharp [2(]]. This 
approach has been very influential and has applications to the IBM in which the same U(5) and SO (5) D SO (3) groups 
also appear. Its limitation is that it is impractical for the description of collective model states of large deformation 
which converge extremely slowly in a spherical U(5) basis. 



III. AN ALTERNATIVE BASIS FOR THE COLLECTIVE MODEL 



Moving between the above solvable limits of the collective model is complicated by the 'rigidity' of two of the limits. 
Both the rigid-rotor model and the beta-rigid gamma-soft (Wilets-Jean) model have delta function components to 
their wave functions and, as a consequence, they are not realizable in the Hilbert space of the vibrator limit except 
as limits of sequences of normalizable wave functions. This limitation expresses the fact that rigidly-defined intrinsic 
quadrupolc moments are unphysical and incompatible with quantum mechanics (as well as relativity theory). I 
therefore consider an alternative basis for the diagonalization of collective model Hamiltonians and show that it can 
lead to convergent solutions for a wide range of collective model Hamiltonians. 

An important feature of the collective model is that its coordinates separate into orthogonal subsets in much the 
same way as those of a single particle in three-dimensional space separate into radial and spherical coordinates. 
Associated with this observation is the fact that the Hilbert space of the collective model is a direct sum 

H CM = M v = H^M) ® HS°( 5 > , (6) 

V V 

of Hilbert spaces labeled by a seniority quantum number v, which each carry an irrcp of a direct product group 
SU(l,l)xSO(5), where SU(1,1) is the group of scale transformations of the beta (radial) coordinate, defined by 
1 = Q-Q, and SO(5) is the five-dimensional rotation group. As a result of this separation of variables, SO(5)-invariant 
collective model Hamiltonians can be diagonalized with just the Lie algebra of SU(1,1) as spectrum generating algebra. 
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Thus, for example, the collective model has two other exactly solvable submodels given by simply adding a 1//3 2 term 
to either the five-dimensional harmonic oscillator Hamiltonian pj or to the five-dimensional analog of the Coulomb 
Hamiltonian j2lj ■ Such models are known in three dimensions as the Davidson [2^ and Kratzer models • Algebraic 
solutions of three-dimensional central-force problems, have been studied widely in terms of the algebra of the direct 
product group SU(l,l)xSO(3) (reviewed in ref. |24j; an historic account is given in Wybourne's book [2^| and an 
overview of the basic methods in ref. jl^)- The natural extension of the algebraic treatment to five-dimensional space 
with dynamical group SU(l,l)xSO(5) is straightforward The fact that solvable central force problems remain 

solvable on addition of a 1/ r 2 term to the Hamiltonian follows from the observation that the commutation relations 
of the SU(1,1) spectrum generating algebras are unchanged, by the substitution V 2 — > V 2 + k/r 2 . This result holds 
in higher-dimensional spaces. 

Following the standard methods for three-dimensional central-force problems , the spectrum and wave functions for 
the five-dimensional harmonic oscillator 



are given by 



and 



H = -^W 2 + 1 -B^ 2 (3 2 , (7) 



(2n + X v )hio, n. = 0,1,2,... u = 0,1,2,... (8) 



where y va is an SO(5) spherical harmonic for an SO(5) irrcp of seniority v, ^ is an associated Laguerre 

polynomal, and A„ = v + 5/2, b = \J h/Bcu. 

The above (standard) results are obtained by use of the SU(1,1) Lie algebra spanned by the operators 

ii = -V 2 , Z 2 = f3 2 , £3 = -i(?-V + |). (10) 

Because the commutation relations of these operator are unchanged if Z\ is replaced by 

^ = -V 2 + ||, (11) 

it follows || that the harmonic vibrational limit of the collective model extends immediately to a model with beta- 
vibrational Hamiltonian 



H(M = = - V 2 + h -± + B^/3 2 , (12) 



2B\ /3 2 J 2 

The only change in the results of the harmonic vibrational limit for the energies {E nv } and wave functions {^ n va} is 
that the relationship = v + 5/2 is generalized to 



A, = l + v /W§) 2 + /?o 4 . (13) 

This is because the replacement Z\ — > Z[ in the SU(1,1) Lie algebra results 0] in a modification of the value of the 
SU(1,1) Casimir invariant from the value X V {X V — 2) = v(v + 3) + | to the value 

\ v (\ v -2) = v{v + 3)+ r ±+ffi. (14) 

The five-dimensional Kratzer model |2l| is handled in a similar way as shown in the Appendix. 

With explicit wave functions, it is a simple matter to compute observable properties of the model states. For 
example, the matrix elements of an electric quadrupolc operator of the form 

Q M = Z(3q M (15) 

with qM taking the values 

q M = cos 7 V 2 M (n) + ^|Sin 7 (P 2 , M (ft) +V\ M {tt)) , (16) 
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factor into products of a matrix clement of (3 between associated Lagucrre polynomials and a matrix clement of a 
v = 1 SO (5) spherical harmonic. Thus, the latter matrix elements are given by elementary SO (5) Clebsch-Gordan 
coefficients as shown explicitly below. 

The properties of a sequence of collective models with Hamiltonian H((3q) and with (3q ranging from zero to a large 
value are discussed in the following section. It is seen there that the above model provides a sequence of analytically 
solvable solutions which progress from the harmonic vibrational [HW(5)]U(5) D U(5) limit to the asymptotic Wilets- 
Jcan [HW(5)]U(5) D [R 5 ]SO(5) limit. 

Subsequent sections demonstrate that the basis of eqn. |J3}, which reduces an SU(1, 1) x SO(5) D SO(3) dynamical 
subgroup chain, provides a much more rapidly convergent basis for the expansion of the eigenfunctions of an arbitrary 
collective model than does a basis that reduces the chain U(5) D SO(5) 2> SO(3). For a collective model Hamiltonian 
that is SO(5)-invariant, the eigenfunctions and spectrum are readily computed using just the Lie algebra of SU(1,1) 
as a spectrum generating algebra. Convergence is then optimized by selecting a value of /3q at the minimum of the (3 
potential of interest. 

For a general collective Hamiltonian, in which both the SU(1,1) and SO(5) degrees are active, expressions are 
needed for the SO(5) spherical harmonics. They can be constructed by the methods of Chacon et al. [2£j. However, I 
follow a simpler construction based on the observation that the SO (5) spherical harmonics are an orthonormal basis 
for the Hilbert space £ 2 (S , 4) of square intcgrablc functions on the four-sphere S4 with respect to the usual volume 
element given, in standard (7, 0, ip) coordinates, by 

dv = sin37d7dfi, (17) 

where dil = sin# d9 dtp. A non-orthonormal basis for C 2 {S/i) is simply constructed as shown in sect. IV 51 This 
basis is then Gramm-Schmidt orthogonalizcd to give the SO (5) spherical harmonics. This method will be deetailcd 
and used to compute a subset of SO (5) Clebsch-Gordan coefficients in a following paper j^. Thus, the facility to 
diagonalize any Hamiltonian on the collective model Hilbert space @ and to compute, for example, B(E2) transition 
rates between its eigenstates, is limited only by the computational power available. An illustrative calculation of the 
spectrum and E2 transitions for an axially symmetric rotor-vibrator model is given in sect. IVlH 

IV. RESULTS FOR THE BASIC SO(5)-INVARIANT MODEL 

Figure ^ shows the potential energy function 

V ^ = ^^ + f) (18) 
of the Hamiltonian l|12|) for (3q = and 8 together with the corresponding ground-state j3 wave functions. The ratio 




2 4 6 8 10 

P (units of b = \J h/Bu; ) 

FIG. 1: The potential energy components of the Hamiltonian 11211 for /3q = and 8 together with the corresponding beta wave 
functions multiplied by /3 2 . The potential is given in units of huo. 



5 



of the (3- width to the mean value of /3, given for these wave functions, by b//?o is proportional to l/Po^/ai. Thus, as 
PoV^J becomes large, the wave functions approach the Wilets-Jean rigid-beta limit. 

The energy- level spectrum, given by eqns. (JSJ) and (|13|) is shown as a function of (3 in fig.[2]for a range of values of the 
deformation parameter j3o- Fig.|21shows that as f3ay/uj — > oo (the Wilets-Jean limit) the energy levels asymptotically 
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FIG. 2: Energies of the five-dimensional Davidson model as a function of /3q. 



approach a harmonic oscillator sequence with the states of each asymptotic level comprising an infinite number of 
states. Low energy levels for specific values of @q are shown in more detail together with _B(E2) transition rates in 
figs. El and H 
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FIG. 3: Energy-levels and B(E2) transition rates for the Hamiltonian i 1 12 1 with /3o = and 2. Energy levels are labeled by 
their SU(l,l)xSO(5) quantum numbers (n, v). Energies and _B(E2) transition rates are given in units such that the one-phonon 
L — 2 excited (0,1) state has an excitation energy of 100 units and its _B(E2) transition rate to the ground state is also 100 
units in the U(5) (flo = 0) limit. 

B(E2) transition rates are computed as follows. Let {|nwcr)} denote the basis states corresponding to the wave 
functions of eqn. © with related to v by eqn. (|Ll|) . The B(E2) transition rate between two SO(5) levels is then 
defined in the standard way for the quadrupole operator Q by summing the squared matrix elements of Q over the 
states of the final level and averaging over initial states. Thus 

B(E2;mVi-*nfVf) = -77— rl^/^/^/IQ/A/l^i^^)! 2 
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FIG. 4: Energy-levels and B(E2) transition rates for the Hamiltonian 11211 with fia = 4 and 6 in units defined in the caption 
to fig. 01 For comparison, (c) is the spectrum obtained by a simple adiabatic approximation with < /3 >= 6.07 as described in 
the text. 



= |^yl(»/«/lllOll|n^>| 2 

= |(n^|||Q|||n /W/ )| 2 , (19) 

where d(v) is the dimension of the SO(5) irrep of seniority v and (n/i?/|||Q|||niiJi) is an SO(5)-reduccd matrix element 
of the v — 1 quadrupolc tensor. Since the wave functions for the states (Infer)} are products of /3 wave functions and 
spherical harmonics, the matrix elements of the Qm = Z(3q~M operators of eqn. I|15|) factor and are given by 

(niVi\\\Q\\\nfV f ) = Z{niVi\(3\nfVf){vi\\\q\\\vf), (20) 

with obvious notation. The (3 matrix element is readily evaluated for any value of /3o while the q matrix element is 
independent of f3o and can be evaluated in the U(5) (/?o = 0) limit. 

In the U(5) limit, the Qm operator can be expressed in terms of L = 2 harmonic-oscillator raising and lowering 
operators 

Qm = Z^I±^L. (21) 

It follows that (OvIIIQIHOu) = and 



(0,t> + l|||Q|||0i;> = Zy/{v + l)/2. (22) 



Therefore, since the (3 integral gives (0,u + l|/3|0«) = y/(2v + 5)/2, we infer that 



(v\\M\v + 1) = (-W^-^- (24) 
y 2v + 5 a(v) 



7 



where the phase factor (—1)^ depends on a choice of phase convention. Thus, for arbitrary Po, 

z 2 

B(E2;n'v - n,v - 1) = ^-^(n'z;|/3|n, « - l) 2 , (25) 

B(E2; n't- - n, « + 1) = ^±^±ll( n 'v\P\n, v + l) 2 . (26) 

(2v + o)d(v) 

It is also follows from the identity (0w|||(5|||0u) = that all quadrupole moments are zero in this model as expected 
for a gamma-soft model. For the present calculations, Z 2 was given the value Z 2 = 200 so that B(E2; 2\ — ► 0i) = 100 
in the U(5) (x = 0) limit. 

It is seen from the figures that, as the value of Po is increased, excited beta-vibrational bands separate and increase 
in energy. In the Wilets-Jcan limit, in which Po^/uj — > oo, the ratio of the beta-vibrational energy to the lowest 
rotational excitation energy diverges. As this limit is approached, a simple adiabatic approximation gives increasingly 
good approximate solution to the results of the Hamiltonian (|12|l as illustrated by comparison of figs. 4(b) and 4(c). 

The adiabatic approximation follows by expanding the expression for the energy JHJl with A„ given by cqn. (|13(1 . for 
large values of po, 



1 



2n + t^2 v ( v + 3) 



huj + I I + const., (27) 



and neglecting terms of order 1/ Po and higher. This approximation corresponds to neglecting the centrifugal coupling 
between the rotational and beta-vibrational degrees of freedom. Perturbation theory shows that this approximation 
is good so long as the rotational energy v(v + 3)/2/?g, for the gamma-soft rotor, is small in comparison with the 
beta-vibrational energy Hlu. 

The effective decoupling of the beta-vibrational and rotational degrees of freedom in the adiabatic limit, is evidenced 
by the independence of the beta wave function on the seniority quantum number v in this limit. A comparison is 
made between the v = 0, 5, and 10 beta wave functions in fig. [S] for /3 = 4, 6, and 8. It is seen that centrifugal 
stretching is substantial for Po = 4 but negligible for Po ~ 8 and v < 10 




FIG. 5: Comparison of the beta components of the model wave functions for Po — 4, 6, 8 and for v = 0, 5, and 10. It is seen 
that the wave functions are shifted to higher values of /3 with increasing v due to centrifugal stretching for small values of Po 
but that as (3oy/v — > oo, the adiabatic limit is approached in which this shifting becomes more and more negligible for finite 
values of v. 



V. MORE GENERAL COLLECTIVE MODEL HAMILTONIANS 



An extension of the above model to admit gamma as well as deformcd-beta equilibrium shapes, is given, for example, 
by adding a gamma-dependent potential to the Hamiltonian i|12[) ; 

6 = 5 ( ~ y2 + $) + l Buj2p2 + v ■ (28) 

To illustrate, I consider the potential 



V(t) = ~x cos 37 



(29) 
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which has an axially symmetric minimum. This potential produces a very tractable model because, to within a 
constant, cos 37 is the S0(5) v = 3,L = spherical harmonic. An even simpler model is obtained by assuming a 
value of Pq^/lo for which the (3 and SO (5) degrees of freedom are adiabatically decoupled. The latter choice is not 
necessary but simplifies the calculations considerably. 



A. The adiabatic map 

For each value of the /3 quantum number n, the SO (5) quantum number v runs over the complete set of nonnegative 

integers v = 0, 1, 2, Moreover, just as the set of SO(3) spherical harmonics {Ylm} spans the Hilbcrt space C 2 (S2) 

of square integrable functions on the two-sphere, so the set of SO(5) spherical harmonics {3^<tm} spans the Hilbert 
space C 2 {Si) of square- integrable functions on the four-sphere (isomorphic to the factor space SO(4)\SO(5)). It 
follows that 



£2 ( S 4)= H SO(5)^ (30) 



v=0 

Observe also that the SU(1,1) Hilbert ' featuring in cqn. © are all isomorphic to the common Hilbert 

space H su ( 1,:l ) = C 2 (R + ) of functions of j3 £ R + (the positive radial line) that are square integrable with respect to 
the measure /3 4 d/3. Thus, we can define an adiabatic map from H^ U( ' 1 ' 1 ^ to the isomorphic Hilbert space Hq 1 in 
which the basis wave functions map 

4>nv l-> 4>n0 = 4>n ■ (31) 

Under this map 

H CM^ H SU(l,l) 0£ 2 ( 3 4) ^ (32) 

The corresponding map of the Hamiltonian is given by the observation that, for (3q large compared to (v + 3/2) 2 , 

K = l + J{v + l? + fit Ao + i«(« + 3). (33) 

Since v(v + 3) is an eigenvalue of the SO(5) Casimir operator C so ( 5 ) ; it follows that the adiabatic map sends the 
Hamiltonian H(/3o) to 

H((3 )^H ((3 ) + ^C so ^, (34) 

where Hq(Pq) is the SU(1,1) Hamiltonian H(/3q) restricted to the scalar v = representation on Hq 11 ^ 1 ' 1 ' 1 . Thus, 
the Hamiltonian Hq(Po) gives a harmonic sequence of /^-vibrational states with energies {nhiu,n = 0, . . . , oo}. The 
Hamiltonian H of eqn. I|28|) likewise maps, in the adiabatic limit, to 

H(/3 ) - i? ad (/3 ) = JT (A>) + ^ C SO(5) + V. (35) 

I consider the Hamiltonian 

H(x) =C so & - X cos3 7 (36) 

as a component of the adiabatic Hamiltonian 

Then, if the energy levels of 7i(x) are given by {£ l}, the energy levels of -ff a d(A)j x) are given by 

Wo 



HM, X ) = H (Po) + 7^2 \c SOi5) - Xcos3 7 l . (37) 



E na L = ntiU + —J £ a L- (38) 
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B. Basis wave functions and matrix elements 

The following construction makes substantial use of the methods of Chacon, Moshinsky, and Sharp |2£} but is 
conceptually and computationally simpler. 

The primary observation is that the Hilbcrt space of square intcgrable functions on the four-sphere is spanned by 
polynomials of the elementary v — 1, L = 2 quadrupolc moment functions: 

<Zm(7,0) = coB72?g fW (n) + ^ sin 7 [^M^ + ^-2,m(^)] • (39) 

To construct a basis, I start by forming a minimal set of angular-momcntum-coupled wave functions, with angular 
momentum projection M = L, from which a complete set of M = L wave functions can be generated by taking 
products of the wave functions in this minimal set. A huge advantage is gained by proceeding in this way because all 
the wave functions generated have good angular momentum and form a complete set of highest- weight wave functions 
for the SO (3) irrcps in C 2 {S A ). 

A suitable minimal set of wave functions is found by examination of the angular momentum content of the L 2 {SjC) = 

(§) v Hi S0 ^ 5 ' ) ' ) space. Knowing the SO(5)— >SO(3) branching rules for the irreps appearing in the collective model |28[ . 
it is possible to arrange the angular-momentum states of C 2 (S A ) into K bands, each having the same sequence of 
angular-momentum states as those of axially-symmetric rotor bands of the same K, and with the sequence of K 
bands being in one-to-one correspondence with those of a sequence of gamma- vibrational bands. More precisely, the 
band-heads of the K bands appear with increasing seniority v in the sequence 

v= 0123356789 10 •■■ 

K= 2 4 6 8 10 ••• .... 

2 4 6 -.. ( 4 °) 
2 4 ••• 

This band structure, albeit with different energies, is precisely that of the axially-symmetric rotor-gamma-vibrator 
model (i.e., no beta-vibrational bands) as can be seen in the spectrum of the Hamiltonian 7i(x) for x — of fig. El It 
is now seen that a complete set of M = L coupled polynomials in is generated by taking multiple products of the 
four generating functions 

$0022 oc <72, (41) 
$0222 oc [q® q] 22 , (42) 
$1000 oc [q<g>q®q]o, (43) 
$0233 oc [q® q ®q}33 • (44) 

Thus, these functions generate a linearly-independent basis of M = L polynomials 

®tKLL = [$0022] ,ll [$0222]" 2 [$1000]*[$0233]" 3 , (45) 

with 

K = 2n 2 + 2n 3 , L = 2n 2 + 2n 2 + 3n 3 , n 3 = orl. (46) 

If a state belonging to a K band is labeled by the index t = i — 1 when it appears in the i'th occurence of K (e.g., 
t = for states of the lowest K = band and t = 1 for states of the first excited K = band), then the polynomials 
of eqn. H45|) are in one-to-one correspondence with the states labeled by the if-band system given above; t is then the 
number of $1000 zero-coupled triplets as in the basis construction of |2Cj . 

The polynomials {$(kla/} do not form an orthonormal basis. However, their overlaps are readily evaluated using 
the inner product for C 2 {S A ) defined by the standard volume element (|17|l . Moreover, since a polynomial QtKLL 
of degree N in q does not contain admixtures of wave functions of seniority greater than N , it is a simple matter 
to sequentially Gramm-Schmidt orthogonalize these polynomial functions to obtain an orthonormal basis of wave 
functions {^ v klm} of good SO(5) seniority v. (Note that the label K of the orthonormal basis is not a good 
quantum number; like the Vergados K label of the SU(3) model it is just a convenient label which makes a useful 
correspondence with the rotor model.) The procedure will be described in detail in ref. p^ . 

In manipulating /^(S^) wave functions, it is convenient to expand them in the form 



even 



(7) V i6 Ja +m [v « m{q) + ^ Lv -^m] ■ (47) 



t>0 
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Thus, a wave function ^ £ £ 2 (S < 4) is defined by a set of functions {^(7), -^2(7), ■ ■ ■}■ In particular, the above 
generating functions are given by 

$0022 ~ {cos 7, sin 7}, (48) 

$0222 ~ {cos 27, sin 27}, (49) 

*iooo ~ {cos 37}, (50) 

$0233 ^ {0, sin 3 7 }. (51) 

The overlaps of wave functions are given in this representation by 

even „ 

= 8ll>6mm> £ / F* Lk {i)F'lM sin3 7 d 7 (52) 

k>0 J 

and matrix elements of the interaction cos 37 by 

even „ 

(* LM |cos3 7 |*^ M ,) = 5 LL ,S M M' / F *lM F 'lM cos3 7 sin37d 7 (53) 

In the present calculation, the needed matrix elements and overlap integrals were evaluated analytically using the 
interactive mathematics computer program 'Maple'. 



VI. E2 TRANSITION RATES AND QUADRUPOLE MOMENTS 

When the adiabatic approximation is valid, the matrix elements of the quadrupole operators Qm = Z(3qM factor 

(VlWQWVl') =Z(f3)(y L \\q\\y L ,). (54) 

Moreover the expectation value (J3), which varies continously with /3 can itself be treated as a free parameter. Thus, it 
only remains to compute the reduced matrix elements of q. This is straightforward using the wave functions contructed 
according to the methods of the previous section. For example, in the SO(5) (x = 0) limit, the v — 0, L = and 
v = 1, L = 2 states of the orthornormal set {\vKLM)} have wave functions 



10000) ~\/t^2 (55) 



|102M) r 

Thus, with qu given by eqn. 139|) . 
and 

it follows that 



15 
16^2 



16tt 2 

cos 7 ^ A/ (") + -j= sm 7 (V 2 2M (Q) + V 2 _ 2M (Q)) 



(56) 



(102||g||000) = 1 (57) 



B(E2:L i ^L f ) = Z*(f3) ^ L ^ff , (58) 



B(E2 : 102 ^ 000) = iz 2 (/3) 2 . (59) 
5 

The calculated results can be used with arbitrary values of the factor Z((3). The transition rates shown in the 
figures are in units such that, in the SO(5) limit (in which x = and v is good quantum number; cf. fig-EJ) 

B(E2 : 102 000) = 100. (60) 

B(E2) values in these units are given explicitly by the expressions 

r(F2-t Winn B{E2 : u - L f ] - ^ L MW L f)\ 2 (M) 

B(E2 . U - Lf) = 100 B(E2:102 ^ 000) - 500 ^-^ . (61) 
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TABLE I: Values of the B(E2 : vLi — ► v — l,L/)/B(E2 : 12i — > 00i) ratio computed exactly as rational fractions. 
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Quadrupole moments are given in the corresponding units by 



Q(L) = V500 (LL, 2Q\LL) ^B= , (62) 

where (LL, 20|LL) is a standard SO(3) Clcbsch-Gordan coefficient. 
In the SO(5) limit, when x = and v is a good quantum number, 

{V+ )^ V ; Lf) = 12||« + 1, L,) (t, + l|||4|||t,) (63) 

where (vLf, 12||u + l, Lj) is a reduced SO (5) Clcbsch-Gordan coefficient and + is the SO(5)-reduced matrix 

element given by eqn. Q24|). Thus, in the SO(5) limit, 

B(E2:v + l,L i ^vL f ) _5(v + l). T 10| , , , n r , 2 



B(E2 : 102 -> 000) (2w + 5) 



wL^^Hw+l,^) 2 . (64) 



Low- lying values of these ratios are given in Tabic [I] Quadrupole moments vanish in the SO (5) limit because q is a 
V = 1 tensor operator and, hence, its matrix elements satisfy the well-known Av = ±1 selection rule. 



VII. RESULTS 



The low-energy spectrum of the Hamiltonian 1|36|) is shown for x = 0, 25 and 50 in figs. and|$l respectively. 
The results were obtained by diagonalization of the Hamiltonian in a basis of 12 states for each angular momentum; 
it was ascertained that this number was sufficient to obtain results of the desired accuracy for x < 50. It can be seen 
from these figures that the spectrum and E2 transition rates progress with increasing x from those of the Wilets-Jean 
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gamma-soft model (at \ = 0) to those of the adiabatic axially-symmetric rotor-vibrator model. In particular, the 
K bands acquire the characteristics of an axially-symmetric rotor model and sequences of excited gamma-vibrational 
bands emerge. This is most evident in the sequence of K = 2, 4 and 6 bands, which appear as one-, two-, and 
three-phonon gamma-vibrational bands, and the move of the first excited K = band-head towards the energy of the 
two-phonon K = 4 gamma-vibrational energy. 




143 143 



— i — I 
100 

L -* — 

FIG. 6: The low-energy spectrum of the Hamiltonian 1361 in the SO(5) limit in which seniority v is a good quantum number. 
Some pairs of degenerate energy levels are separated slightly for clarity. The B(E2) transition rates shown are in units such 
that B(E2; 12i — * 00i) = 100. A more complete set of transition rates is given (precisely) in table [I] Quadrupole moments 
vanish for states of good seniority. 

It is instructive to compare the results of fig. |S1 with the sequence of rotational bands expected in the standard 
axially symmetric rotor model shown in fig. [5] 

The comparison shows that the rotational dynamics are adiabatic relative to the gamma-vibrations for the lowest 
energy states but that significant departures become evident at higher energies. Increasing the value of x would surely 
reduce the rotational-vibrational coupling and give results closer to those of the adiabatic rotor model. However, by 
the same token, the gamma-vibrational bands would rise to higher energies. It is of interest to see if some of the effects 
of rotation- vibration coupling, e.g. the non-rotational behaviour of some of the E2 transitions in the two-phonon K = 4 
gamma band, arc realized in rotational nuclei having low-energy gamma bands. It would also be interesting to see if 
the results of fig. |S] can be reproduced more accurately in the adiabatic rotor model by taking the rotation- vibration 
coupling into account in perturbation theory. 

VIII. SUMMARY 

The Wilets-Jean model is an algebraic model with an [R 5 ]SO(5) D SO(5) D SO(3) dynamical subgroup chain. 
However, as observed above, the states that reduce this dynamical chain are strictly zero-width beta-rigid states. 
Complete rigidity is unphysical. Nevertheless, it is possible for the WJ model to work well even when the beta 
fluctuations of physical states are relatively large. This is because, it is hard to measure the beta-width of physical 
states; indeed it is invariably unclear as to the nature of excited K = excitations. Many intrinsic excitations of a 
nucleus can give rise to excited K = bands and even the collective model has two-phonon K = gamma bands. 

To avoid the problems associated with the beta-rigidity of the WJ model, we considered in this paper a more physical 
collective model which, in its gamma-soft limit, has an SU(1, 1) x SO(5) D SO(5) Z) SO(3) dynamical subgroup chain. 
For large values of the parameter f3oy/uj : the low-energy states of this model acquire properties that are essentially 
indistinguishable from those of the less physical WJ model. In the WJ limit of this model, the beta vibrational 
frequency is infinite for a finite value of (3q and, hence, beta vibrational excitations are not observable. However, 
neither an infinite beta-vibrational frequency nor a zero beta-width for collective wave function is necessary to obtain 
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the results of the WJ model. All that is needed is an adiabatic decoupling of the rotational and beta-vibrational 
degrees of freedom. The signature of such decoupling is that the beta wave functions in the SU(1, 1) x SO(5) model 
become essentially independent of the excitation energy for the range of energies of interest. That such is the case 
when the rotational energies arc small compared to the one-phonon beta-vibrational energy is illustrated in fig. [S] 

Relaxing the rigidity of the WJ model and replacing its [R 5 ]SO(5) dynamical group with SU(1, 1) x SO(5), results 
in a more physical collective model. More importantly, it makes it possible to construct a basis in which arbitrary 
collective model Hamiltonians can be expanded in a rapidly convergent manner. For example, the solutions of the 
Hamiltonian 13 7|) given in figs. KilSI were obtained by diagonalization of 12 x 12 matrices. For values of \ > 50, 
the dimensions would have to be increased to ensure accurate results. But clearly, much larger dimensions can be 
accommodated with available computers. 



APPENDIX A: THE FIVE-DIMENSIONAL KRATZER MODEL 

The five-dimensional Kratzer model, proposed by Fortunato and Vitturi pH , has a similar algebraic solution to 
that of the Elliott-Evans-Park model [j] . A straightforward extension of the methods used for the hydrogen atom to 
the Hamiltonian 

H = V 2 , (Al) 

of a five-dimensional Coulomb problem gives the energy levels 



k 2 B 

and wave functions 
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FIG. 8: The low-energy spectrum and B(E2) transition rates for the Hamiltonian 1361 with \ = 50. The units are as defined 
in the caption to fig. For \ non-zero, seniority is no longer a good quantum number and quadrupole moments are no longer 
zero. 

where now A r = 2t + 4 and b = H 2 /(2kB). 

These results are obtained by use of the SU(1,1) Lie algebra spanned by the operators 



Zi = -2/3V 2 , Z 2 = 2/3, Z 3 = -2i(q ■ V + |) . 
Again it is seen that, the commutation relations of this Lie algebra are unchanged if Z\ is replaced by 



(A4) 



(A5) 



It follows that the above results for the Kepler Hamiltonian, extend immediately to the five-dimensional Kratzer 
model with Hamiltonian 



H(e) 



2B 



V 2 



pJ (3 



(A6) 



However, now the relationship A r = 2r + 4 is replaced by 



A r = l + 2J(T + l) 2 +ffi. 



(A7) 



This follows because the replacement Z\ — > Z' x in the SU(1,1) Lie algebra modifies the value of the SU(1,1) Casimir 
invariant from the value A T (A T — 2) = 4t(t + 3) + 8 to the value 



A T (A r - 2) = 4r(r + 3) + 8 + 4/3^ 



(A8) 
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FIG. 9: The low-energy spectrum and B(E2) transition rates given by the axially symmetric rotor model with band-head 
energies and moments of inertia adjusted so that the dashed energy levels are fitted to those of fig. [S] The intrinsic quadrupole 
moments of the model were also adjusted to fit the E2 transition rates shown in smaller type. 
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